Scaling of nascent nodes in extended s-wave superconductors 
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We analyze the low-energy properties of superconductors near the onset of accidental nodes, i.e. 
zeroes of the gap function not enforced by symmetry. The existence of such nodes has been motivated 
by recent experiments suggesting a transition between nodeless and nodal superconductivity in iron- 
based compounds. We find that the low-temperature behavior of the penetration depth, the specific 
heat, and the NMR-NQR spin-lattice relaxation rate are determined by the scaling properties of 
a quantum critical point associated with the nascent nodes. Although the power-law exponents 
are insensitive to weak short-range electronic interactions, they can be significantly altered by the 
curvature of the Fermi surface or by the three-dimensional character of the gap. Consequently, the 
behavior of macroscopic quantities near the onset of nodes can be used as a criterion to determine 
the nodal structure of the gap function. 



After three years of their discovery, the symmetry of 
the superconducting state of the iron-based compounds 
is still under intense debate Although many early 
probes suggested a fully gapped state, recent experi- 
ments have found strong evidence for the presence of gap- 
less quasi-particle excitations, suggesting the presence 
of nodes in at least some regions of the phase diagram 
(3-[lo|. Interestingly, in overdoped Ba (Fei_2;rMa;)2 As2 
materials, with TM = Co, Ni, Pd, penetration depth 
[2, 3], specific heat Q, and thermal conductivity mea- 
surements [H, 01 indicate a transition from a nodeless to 
an anisotropic nodal state, that does not seem to be re- 
lated to changes in the Fermi surface (FS) topology. In 
underdoped {Qa-i-x^x) Fe2As2 compounds, where super- 
conductivity (SC) may coexist with the magnetically or- 
dered phase, recent thermal conductivity and penetration 
depth measurements also show a doping-induced onset of 
nodal quasi-particles l^, 11|- 



One possible scenario to explain these nodeless-to- 
nodal transitions is a change in symmetry from a SC state 
with finite gap everywhere on the FS, e.g. an isotropic 
s-wave pairing state, to a state with symmetry-enforced 
nodes, such as a d-wave state. In this case, a line of 
finite temperature phase transitions would separate the 
nodal and nodeless SC states. Indeed, for a range of 
band structure parameters, some theoretical models find 
nearly degenerate s-wave and d-wave solutions of the gap 
equations associated with purely electronic pairing inter- 
actions [12 , 13J. Another scenario, not associated with 
a finite-temperature phase transition, is the emergence 
of accidental nodes, i.e. zeroes of the gap function at 
positions of the FS unrelated to a particular represen- 
tation of the tetragonal symmetry group. In the iron 
pnictides, one frequently discussed gap function A (k) is 
the extended s-wave, with different signs along different 
FS sheets The location of its zeroes in momen- 

tum space is determined by energetic arguments rather 
than symmetry requirements. Therefore, changing the 
size of the FS sheets or modifying quantitative details 
of the pairing interaction can cause such a zero to touch 
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Figure 1: Accidental nodes in spherical (a) and cylindrical (b) 
Fermi surfaces. The point and the full line denote the 5 = 
nascent node in each case, whereas the dashed lines denote 
the line nodes that develop for 5 < 0. d\\ is the number of 
directions on the FS with curvature. In (c) we present the 
T-5 phase diagram with the crossover regime (shaded yellow 
region between the dashed lines) and the scaling properties of 
the variation of the penetration depth AA. 



and eventually cross one of the FS pockets, leading to a 
nodal state [Tsl-^. 

At finite temperatures, the onset of accidental nodes 
does not imply a phase transition. However, it is charac- 
terized by a quantum critical point (QCP) that separates 
a regime with finite-ranged quasi-particle excitations, in 
the fully gapped state, from a regime with algebraically 
decaying excitations, in the nodal state. Thus, general 
crossover arguments of QCP should be useful to charac- 
terize the behavior of macroscopic quantities near the on- 
set of accidental nodes. Consider, for example, the varia- 
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tion of the penetration depth AA = A (T) — A (0) (we only 
consider the clean case here) : in the gapped state one ex- 
pects activated behavior, AA cx exp (— A/T), whereas in 
the case of line nodes holds AA = AT, with constant 
A. Naively, one expects that the transition between the 
two regimes is characterized by a power-law AA = A'T" 
with exponent 1 < n < oo intermediate between the two 
regimes. For instance, if the accidental nodes first touch 
the FS at a single point, one would expect to observe 
n = 2, characteristic of systems with symmetry-enforced 
point nodes (see Ref.|29|] for a review). 

In this Brief Report we analyze the crossover behavior 
near the onset of accidental nodes and how it depends 
on the shape of the Fermi surface. Among others, our 
results demonstrate that even if nascent nodes originate 
around a single point, they behave qualitatively different 
from symmetry-enforced point nodes. Our main results 
are summarized in Fig.l: in the case of a spherical FS, 
at the onset of accidental nodes, we obtain a linear-in-T 
penetration depth (n = 1), but with a coefficient that is 
precisely half of the value that one finds for line nodes, 
A' = A/2. Consequently, after the onset of accidental 
nodes, AA displays a crossover between two linear-in- 
T regimes with different slopes. Similar behaviors hold 
for other quantities, such as the specific heat C and the 
NMR-NQR spin-lattice relaxation rate Ti. In the case 
of a perfectly cylindrical Fermi surface we find, in agree- 
ment with earlier calculations [s^, that n — 1/2, i.e. a 
behavior more singular even than for line nodes. De- 
noting by the number of directions on the FS with 
curvature (dy = 2 for a spherical FS and = 1 for a 
perfectly cylindrical FS), our results for the scaling be- 
havior are expressed as AA cx T''"/^, C oc T^+'^i\/^, and 
T^^ cx T^^'^ii . In order for the more exotic regime of 
cylindrical FS (cZy = 1) to be observed, the energy asso- 
ciated with the FS curvature must be negligible on the 
scale of the SC gap Aq ~ 10 meV, which is a very restric- 
tive condition. Yet, this sensitivity with respect to the 
shape of the FS might enable one to distinguish between 
distinct nodal structures of the SC gap. In what follows, 
we will derive these results and show that they are robust 
against weak short range electron-electron interactions. 

To describe the onset of accidental nodes we consider 
the angle-dependent gap function on the FS: 



and TT. Our results do not depend on the specific form 
of A {9) given in Eq[l] rather, they apply generally to 
systems where, close to the nodes, the gap can be written 
as A {0} ~ Ao{9 + 6*0+) (6* -9o-)- In this language, the 
QCP is characterized by the merging of the two nodes, 
giving rise to a quadratic dispersion A {9) ~ Ao {9 — 9o)^. 
In the context of the iron pnictides, 9 in EqlT] could be 
an angle relative to the center of the electron pocket. If 
one considers instead accidental nodes in the hole pocket, 
one simply replaces cos (20) by cos (40), which does not 
change the results obtained here. 

The low temperature behavior of numerous observables 
2^ . 31 1 is determined by the density of states of Bogoli- 
ubov quasi-particles: 



N (w) = iVt 
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where Nq is the normal-state density of states at the 
Fermi level and cKl refers to the angular integration. 
With N {lo), it is straightforward to obtain the tempera- 
ture dependence of the penetration depth, the specific 
heat, and the NMR-NQR spin- lattice relaxation rate. 
Denoting the Fermi function by / [uj) — {e^^ + l) , 
with (3 — {ksT) ^ , one obtains: 
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C cx 



duN (w) 



duN {uj) (3uj^ 
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Band structure calculations as well as ARPES mea- 
surements reveal that the FS of the iron pnictides are 
anisotropic, similar in shape to warped cylinders (s^ . 
Although small, this curvature of the FS becomes very 
relevant for the low energy properties of interest here. 
Thus, close to the accidental nodes, the FS is similar 
in shape to an ellipsoid which, upon rescaling, can be 
described by a sphere, yielding d^l = (47r)~^ sm9d9d(p in 
Eq. [2l In this case the accidental line node becomes ellip- 
tical in shape, starting at an isolated point (see Fig. la). 
The resulting density of states is given by: 



A(6>) _ 1 + 6 1-S 
Aq ^ ~ 2~ 



cos (26') 



(1) 



Here Aq > is the maximum amplitude of the gap, 9 
is the angle measured relative to the center of the FS, and 
S is the tuning parameter. For S > the superconductor 
is fully gapped with minimum gap amplitude AqS. For 
(5 < 0, the gap function changes sign at nodes located at 
angles not fixed by symmetry 9o± — ±i arccos ^y^l^ , 
and its most negative value is AgfJ < . The special 
case S = corresponds to the onset of nodes at = 
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In Fig. 2, we show N (w) for 6 larger, smaller and equal 
to zero. At low energies |a;| <C Aq and near the onset of 
nodes (i.e. for \5\ ^ 1), we obtain: 
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Figure 2: Density of states A'^ as function of energy u in the 
case of a spherical Fermi surface for different values of S. The 
insets show the dependence of the gap A on the polar angle 



where g (x) = 1 for 
for 



iarccotV — 1 for \x\ > 1. 
point nodes at S — behave differently than the symme- 
try enforced point nodes occurring in the polar state of 
triplet superconductors (29l . 31 1. In the latter, the gap 



< 1 and g (w) = ^ + 
Notice that the nascent 



magnitude near the node varies as |0 — ^qL ^ opposed 
to {9 — 6*0)^ discussed here, giving rise to N (w) oc in- 
stead of the N (cj) cx oj found here. A similar result was 
discussed in Ref. [33], in the context of hybrid pairing 
states of uniaxial superconductors. 

Because of the shallow rise of the gap near the nascent 
nodes, the slope of the density of states for 5 = is 
only half of the value for 5 < 0, where the gap actually 
changes sign. Consequently, for (5 < a crossover takes 
place between two regimes that are linear in uj, with a 
low-energy slope twice as large as the higher-energy slope. 
In the fully gapped regime, d > 0, the density of states 

-<5 ' 



vanishes as Nn 



near the threshold, but 



1-.5 V Ao 

soon recovers the linear critical behavior at higher en- 
ergies. This r = crossover regime is also manifested 
at finite temperatures, as we show in Fig. 3, where the 
T-dependence of the penetration depth is displayed for 
different values of the tuning parameter 6. To probe the 
low-energy behavior of N (ui), low temperatures are re- 
quired. Therefore, the observation of this crossover, or 
of a similar crossover in the slope of C/T or 1/T^Ti, 
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Figure 3: Variation of the penetration depth AA = A (T) — 
A (0) as function of temperature T for several values of S. 
Notice in all curves the crossover at higher temperatures to 
the Unear-in-T dependence associated with the nascent nodes 
(solid red line). At low temperatures, AA decays exponen- 
tially for the fully gapped case (solid blue line) and linearly 
for the nodal case, but with a different slope (dashed green 
line). Here we used the values S = 0.025, 0.01, 0, -0.005, 
-0.01, -0.015, -0.02, and -0.025. 



would be a strong evidence for accidental nodes and thus 
extended s-wave pairing. 

For completeness, we also discuss the regime of per- 
fectly cylindrical FS. In this case the angular integration 
becomes dfl — {2tt)^^ dd and two nearly parallel line 
nodes emerge from a single line of nascent nodes with 
A {9, kz) ~ Ao (6* — ^o)^ independent of (see Fig. lb). 
The low-energy density of states is then given by 
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Nn 
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with g (x) 



-3/2, 



1 

1-1/2 



for 
for 



1, 



9ix) 

-1 



1 < < \6\ , and 



logarithmically diverging g (x) = \/2|(5| /tt^ log (j^pi 
close to x = 1. The surprising aspect of this result is 
that the power in N (uj) at the QCP is lower even than 
the case of line nodes, a result that was obtained earlier 
in Ref. (s^l- Thus, right at the transition to a fully 
gapped regime there are more low energy quasi-particle 
excitations compared to the nodal regime, again a 
consequence of the shallow rise of the gap for arbitrary 
kz- Since the actual FS is not perfectly cylindrical, 
one expects a crossover to the linear in uj behavior for 
a spherical FS, Eq. [S] This crossover will depend on 
details of the fc^-dispersions of the FS (as observed in 



m) 



|32l |) and of the SC gap (as observed in [34 
particular, on how the curvature of the FS changes as 
function of energy when compared to the changes in 
the SC gap. In Ref. Reid et al. suggest the two 
nodal structures associated with Eqs. [S] and IH] as the 
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main candidates to explain their thermal conductivity 
measurements in overdoped Ba (Fei_2;Coa;)2 As2. The 
different crossover behaviors of several macroscopic 
quantities that results from Eqs. [S] and [HI provide an 
important criterion to decide between the two scenarios. 

The results above can also be obtained from scaling 
arguments of the low-energy action associated with the 
QCP. The free fermion action for ^ = is 



So = / v4 



; — vkj_ 



where we used the Nambu spinors ipk = ( c^^ 




-ki 



(7) 



to 



denote the nodal quasi-particles. Here k — (fcj^, k|| , iw) 
stands for momenta perpendicular and parallel to the 
FS measured from the location of the nascent node, 
with J^, = TJ2„ J dk±d'^«k^^. Upon scaling k± bk±, 
— )• and u) — >■ buj, the action becomes scale in- 

variant if one uses the rescaled field ipk = b'^4''k' with 
K = (6 + c?||) /4. Thus, we find the scaling dimension of 
short range electron-electron interactions to be — (i||/4. 
Since this scaling dimension is negative, weak short range 
electronic interactions are irrelevant in the sense of the 
renormalization group theory and should not lead to an 
instability. It remains to be analyzed the effects of long 
range interactions, such as quantum phase fiuctuations. 

Since the scaling dimension of the tuning parameter 5 
is unity, we obtain the low-energy scaling form for the 
density of states: 



N{u:,5) = h''^\\/^N {hu:,h5). 



(8) 



For (5 = 0, we can fix the arbitrary scaling parameter 
h = ijj'^ which yields N (w, 0) oc uj'^wl'^. For finite 5 < 0, 
it holds that the slope A {5) = dN (ui, S) /du!\^^Q behaves 
as A{S) oc These results are fully consistent 

with the explicit analysis presented above. 

In summary, we analyzed the scaling behavior near the 
onset of accidental nodes in extended s-wave supercon- 
ductors, and how it depends on the curvature of the FS. 
The low-T properties of the system can be analyzed in 
a simple scaling theory that is governed by a QCP sepa- 
rating a nodal and a nodeless region, without symmetry 
breaking. A somewhat similar behavior was proposed 
in extreme type-II superconductors, where the spectrum 
of Landau levels becomes gapless beyond a threshold 
magnetic field [36]. Our theory predicts that the low- 
temperature slope of the linear quasi-particle density of 
states at the onset of accidental nodes is reduced by a 
factor of 1/2 compared to the regime with well estab- 
lished nodes. This leads to an experimentally observable 
crossover in a number of physical observables, allowing 
one to distinguish whether nodes are accidental or sym- 
metry enforced. 
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